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A MATHEMATICAL MODEL OF COPPER CORROSION 

F. CLARELLI, B. DE FILIPPO, AND R. NATALINI 

Abstract. A new partial differential model for monitoring and detecting 
copper corrosion products (mainly brochantite and cuprite) is proposed to 
provide predictive tools suitable for describing the evolution of damage 
induced on bronze specimens by sulfur dioxide (SO2) pollution. This 
model is characterized by the movement of a double free boundary. 
^ 1 Numerical simulations show a nice agreement with experimental result. 
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Deterioration of copper and bronze artifacts is one of the main concerns 
for people working in cultural heritage [25]. More specifically, bronze, 
a copper-tin alloy, has been widely employed for daily-use and artistic 
purposes from the Bronze Age up to present. Conservation studies are 
based on a knowledge of the environmental conditions to which copper 
and copper alloys may be exposed and include all the information on the 
material technologies and the nature of the corrosion films or patina, which 
cover their surfaces. In particular a significant effort has been devoted to 
study the corrosion due to environmental conditions, such as temperature, 
^ ! moisture, concentration of pollutants It3l [T9l [TTll . Although in recent years 

<3\ " air pollution in European urban areas has decreased considerably, there 

still remain concentrations of pollutants such as sulfur dioxide (SO2) from 
combustion of fossil fuels, being one of the most important factors in the 
deterioration of bronze. Indeed SO2, mixed with water vapor, reacts to 
produce sulfate acid (H2SO4), which causes corrosion phenomena on copper 
surfaces and produces several corrosion products as basic copper sulphates, 
such as antler ite, posnjakite, brochantite llT4ll . The latter is the final product 
r> ! of several reaction steps, which can be approximated by two main chemical 

reactions: cuprite formation, which occurs after a few weeks of exposure to 
atmospheric conditions, and brochantite formation, which the final reaction 
step (for more details on chemical background see Section |2TT) IfTOH . 

The complexity of corrosion processes creates the necessity for 
a quantitative model approach to develop predictive tools, which 
simultaneously provide both quantitative information as well as 
simulations of the various processes involved. These methods, similar 
to those introduced in QQ, are useful for the monitoring and detection 
of surface alterations even before they are visible, making it possible to 
determine optimal intervention strategies. In this paper we introduce a 
new partial differential model, which is used to describe the evolution of 
damage induced on a bronze specimen by atmospheric pollution. It is 
based on fluid dynamical and chemical relations and it is characterized by a 
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double free boundary: one between copper and cuprite, the other between 
cuprite and brochantite. Its calibration has been elaborated according to 
the experimental results in Q. The paper is organized as follows: in the 
Section 12.11 we analyze the main chemical corrosion phenomena and in 
Section 12.21 we review the main mathematical models already proposed 
in literature. ...The section [3] is entirely focused on the description of the 
model's equations and the numerical schemes used, meanwhile in section 
H] we describe the experimental setting and the calibration of the model. 
Finally, in section [5J we present the main results related to the simulations 
produced by our model and in section|6]the related conclusions. 

2. Modeling backgrounds 

2.1. Chemical backgrounds. Most of the oxidation processes occurring 
on bronze metal artifacts, when exposed to environmental conditions, 
are electrochemical and involve interactions between the metal surface, 
the adsorbed moisture and various atmospheric gases (SO2, CO2, NO x , 
hydrocarbons) Il22l . Electrochemical corrosion processes in electrolyte and 
condensed moisture layers have been the subject of extensive studies, based 
on numerous and different approaches ||28l l24l [5]|. When exposed to the 
atmosphere, copper and its alloys form a thin layer of corrosion, from 
brownish-green up to greenish-blue colors, which is designated as patina. 
In the case of copper in a low pollutant levels atmosphere, a native cuprite 
(copper(I) oxide or C112O) film of approximately a few nanometers thick, 
protects the metal surface from further oxidation. The general reaction is 
well described in literature, where, in aerated solution, copper can dissolve 
electrochemically forming copper(I) oxide formation, due to the reaction of 
copper with oxygen. It is represented by the following schematic reaction 
®: 

2Cu + ^0 2 -» Cu 2 0. 

However, in an aggressive environment, like urban atmosphere, the 
protective nature of this oxide layer is altered and there is the formation 
of a non-protective, multi-component, tarnish layer. When the copper 
is exposed to humidity and sulfur dioxide, three kinds of basic copper 
sulfate hydroxide are mainly produced: Brochantite CiiiSO^OH)^, which 
is a well-known patina constituent, or other similar products like Antlerite 
or Posnjakite. In the following we focus our attention to the formation of 
Brochantite, which is the main observed product. Gaseous sulfur dioxide 
and sulfate particles are deposited on the electrolyte on the cuprite. Their 
deposition reduces the pH of the adsorbed water, and this promotes the 
dissolution of cuprous ion (Cu + ) and its oxidation, thus forming cupric 
ions (Cu 2+ ). In detail, copper (I) ions in solution disproportionate to give 
copper(II) ions and a precipitate of copper Q}: 

(1) 2Cu^j - cu\; q) + Cu (s) . 

When the cupric and sulfate ion concentrations in the electrolyte are high 
enough to form brochantite, this phase starts to precipitate on the cuprite 
1H[16]|. In |20ll it is indicated that, in the initial oxidation process, cuprite 
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formation is followed by posnjakite, as a precursor phase to brochantite, 
see also l27l , but we are going to neglect this intermediate transformations, 
due to the elevated speed of the reaction with respect to the time scale of 
the mathematical model which we are going to present. 

2.2. Existing mathematical models. Graedel and his collaborators l|T5ll26ll 
have studied atmospheric copper sulfidation at AT&T Bell Laboratories 
in both experimental investigation and physical or mechanistic model 
development; furthermore a systematic investigation of copper sulfidation 
kinetics has been performed. In the paper 1261 Tidblad and Graedel 
developed a model able to describe the SO2 copper corrosion. Their work is 
based on aqueous chemistry and without considering spatial dimensions. 
In the paper [21], they introduce a model which describes the corrosion 
of copper exposed to moist air with a low SO2 concentration. Here, they 
consider four stages in the development of a corrosion patina: the metal, a 
non-protective oxide film which has a high ionic transport property, an outer 
layer of corrosion products which can permit the penetration of water and 
gases and an external absorbed water layer. More recently, Larson proposed 
in f\7\ a different model that describes the atmospheric sulfidation of copper 
by proposing a physical copper-sulfidation model that includes four distinct 
phases: the substrate metal, a cohesive cuprous sulfide (CU2S) product layer, 
a thin aqueous film adsorbed on the sulphide and the ambient gas. Larson 
postulated that transport through the sulfide layer occurs via diffusion and 
electromigration of copper vacancies and electron holes. Later, in the 1990s, 
copper sulfation by SO2 was investigated by Payer et al., who focused their 
attention on the early stage of corrosion in moist air (75% RH at 25° C) with 
a sulfur dioxide concentration of 0.5%. The techniques employed (SEM, 
AES and TEM) allowed for an analysis and characterization of the oxide 
film (composed of cuprous oxide, copper sulfate and sulfide) on copper 
surfaces and for the mechanism of evolution of the corrosion chemistry 
to be described 13). In the last decade the effect of sulfur dioxide in acid 
rain on copper and bronze has been investigated. Robbiola and his co- 
workers studied the effect of its cyclic action on bronze alloys, underlining 
the different types of patina formed in "sheltered" and "unsheltered" areas 
of bronze monuments EH5J. 

Subsequently, Larson refined the copper-sulfidation model by focusing 
on the transport of charged lattice defects in a growing C112S product layer 
between the ambient gas and the substrate metal. As previously, this 
transport is postulated to occur via both diffusion and electromigration. 

3. The model 

In this section, we aim to introduce a mathematical model able to describe 
the corrosion effects on a copper layer, which is subject to deposition of SO2. 
The present model is based on the mathematical approach used in (6). We 
assume to have a copper sample on which is formed a non protective oxide 
layer (C112O), and, over this layer, a corrosion product (brochantite) grows. 
Over the brochantite layer is assumed to be the atmospheric air with SO2. 
An example of these three layers can be seen in Figured] 
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Cuprite 




Figure 1. Example of cuprite and brochantite deposition on a copper sample. 
The reaction producing cuprite CujO, can be approximated by 



Namely, two moles of copper combined with one-half mole of oxygen 
produce cuprite. 

Laboratory tests, with high concentration of SO2 and high relative 
humidity near to 100%, show that brochantite is the primary product of 
the reaction, thus it can be considered as the final state reached by the 
whole reaction. 

The overall simplified reaction of brochantite formation can be 
approximate by the following reaction ((3J) 



where two moles of cuprite combined with three moles of water and three- 
halves of oxygen produce one mole of brochantite. 

In the following, we assume that these two reactions are instantaneous; 
thus, we obtain a sharp free boundary between cuprite and the unreacted 
copper, due to the reaction (0, and a second free boundary between 
brochantite and cuprite, due to the reaction @). 

Assuming these two reactions as instantaneous, the effective time of 
reaction is implicitly included in the diffusivity coefficients. 

3.1. Swelling. We indicate the copper consumption by a(t). The production 
of C112O on the boundary between copper and the oxide layer proceeds since 
the water and the oxygen diffuse through the oxide layer. On the upper 
boundary of cuprous oxide, a brochantite layer begins to form. It is assumed 
that, in the climatic chamber at a 100% of RH, a thin film of water (with SO2 
dissolved) is formed over the cuprous oxide, and it plays an important role 
in both reactions. By these assumptions, the consumption of CU2O and the 
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production of brochantite on the upper boundary between oxide and water 
film occurs. The volume of cuprite consumption is b(t). 

The transformation of copper into CU2O, such as the transformation of 
cuprite in brochantite, are accompanied by a volume change (swelling rate). 
The swelling rate can be calculated easily because the molar ratio in reaction 
(Q} between Cu and CujO is 2 : 1. Thus two moles of copper change into one 
mole of CuiO, and a different volume of the new matter formed is obtained. 
The swelling of reaction |(T]) is 

(4) d e - -copd; 

where a e is the swelling of reaction (D, and 

represents the expansion volume ratio; \i c and \i v being the molar density 
(moles /cm 3 ) of Cu and C112O respectively. It is assumed that [i c and that \i v are 
constant (i.e. they are homogeneous materials). Under these assumptions, 
if a e (0) = a(0) = 0, then it is possible to conclude that 

(6) a e (t) = -co p a(t), 

and the thickness of cuprite layer h p (t), is proportional to the copper 
consumption a(t) 

(7) h p (t) = {l + co p )a(t). 

On the external boundary of cuprite, we assume that SO2 reacts with 
C112O in presence of water, see reaction ||3). Here, we have that 2 moles of 
cuprite change in 1 mole of brochantite, wasting 1 mole of SO2, 3 moles 
of H2O and 3/2 moles of oxygen. Thus, the brochantite layer grows on 
the external boundary of cuprite, and we indicate the external boundary of 
brochantite with y(t). 

As in equation ©, the swelling rate of brochantite is given by 

(8) b e {t) = -a) b b(t); 
where b e (t) is the swelling of reaction ((3]), 

(9) <o i = £--l, 

and is the constant molar density of brochantite. 

The cuprite consumption b(t) is referred to a system of reference at rest, 
but we have to take into account the moving cuprite boundary. The physical 
boundary between cuprite and brochantite is given by 

(10) p(f) = 6(9 + ^(9, 

where b(t) is the consumption of cuprite (assumed to be positive such as 
a(t)), then a e < 0, see eq. I©. Thus, the variation of the overall external 
boundary y is given by 

(11) y = -a)pA(t) - ojj,b(t). 

Summarizing, the geometry of our problem is one-dimensional, and we 
have 4 regions (e.g. see figure |2]): 
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(1) Copper (inner region). 

(2) Cuprite C112O, between a(t) and j}(t). 

(3) Brochantite Cu^OH^SOi, between /3(f) and y(t). 

(4) Water film and air on the external side of y(t). 



Cuprite 




Thickness (cm) 



Figure 2. Example of cuprite and brochantite theoretical growth (in 
time) on a copper sample. 



3.2. Equation of the model. It is well known that relative humidity plays 
a key role in regulating the speed of oxydation and sulfation. It has been 
observed that when relative humidity exceeds some threshold, then SO2 
reacts completely such as the oxidation of copper happens with full speed. 
We can interpret this phenomenon as follows. According to eq. ((3]), when a 
molecule of SO2 comes in contact with cuprite, it reacts if three molecules 
of H2O are available at the same point (we suppose that there is always 
enough O2). To be more precise, this is true only if vapor condenses on the 
unreacted specimen surface forming a liquid film, and this happens with 
high relative humidity values. 

The brochantite formation (eq. has been assumed to develop on 
the cuprite layer, due to the SO2, H2O and O2, which move through the 
brochantite layer and react with C112O. This reaction implies a wasting of 
C112O, with the formation of a new layer of brochantite. Summarizing, the 
region of brochantite is given by y(t) < x < |3(f) and the region of CU2O is 
j5(f) < x < a. 

Let's denote the concentration of SO2 in the pores of brochantite by S, 
coming from the external air, such as the water concentration indicated by 
W and the oxygen concentration by O throughout brochantite and by G 
throughout cuprite. 

The flow of SO2 relative to air is governed by Fick's law. Thus, in the 
frame of reference where copper is at rest, the SO2 flux has the following 
expression 

(12) /s = n h |-°s^ - Sco p d - Sco b bj = n b |- D s^ + Sfj , 

where the first term refers to SO2 diffusion in the brochantite layer, D s is 
the diffusivity, the second and the third terms on the left side refer to the 
swelling caused by cuprite and brochantite formation respectively. 
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Hence the mass balance of SO2 in the brochantite layer y(t) <x< j}(t) is 

dS „ d 2 S dS n 

(13) *- D '^ + ''^ = 0; 

The value of S at the external boundary y(t) is the environment SO2 
concentration S a (t), which is a known function of time 

(14) s( 7 (t)) = s a (ty, 

We assumed that SO2 reacts totally with C112O at the front j5(t), thus we have 

(15) S(fi(t)) = 0. 

Now, we need of further condition (first free boundary). Since the SO2 flux 
of moles at the boundary j3(£) is proportional to C112O moles consumption, 
we have 

DsdS lp l'U 

(16) - n h — = b; 

v ; b M s dx 2M p 

where M s , M p are the molar weight of SO2 and C112O respectively, p p is the 
mass density of cuprite. 
Similarly, the water flux is 

(17) J w = n b \-D w — - Wco p d - Wa) b b\ = n b \-D w — + Wyl, 

where D w is the water diffusivity Thus, the mass balance in the brochantite 
layer is 

, . dW d 2 W dW 

(18) ^- D »^r + ' > *= a 

The value of W at the external boundary y(t) is the environment water 
concentration W fl (f) 

(19) W(y(t)) = W a (t)} 

Since some moles of water are wasted by the reaction (O, at the front f$(t) 
we have 

(20) l?!L = 3PL b + n b ^b. 



Finally, the oxygen flux is 



dO 



(21) J = n b \-D -^ + Oy\, 

where D is the oxygen diffusivity. Thus, the mass balance in the brochantite 
layer is 

(22) 1F- D ^ + ^- - 

The value of O at the external boundary y(t) is the environment oxygen 
concentration O a (t) 

(23) 0(y(t)) = O a (t). 
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Since the oxygen is also wasted by the reaction ®, at the front |3(£) we have 
(24) A = !^ + n °i. 

3.3. Cuprite layer equations. We assumed that the reaction between 
oxygen and copper occurs on the copper-cuprous oxide boundary a{t). 
Here, we indicate oxygen by G, just to avoid confusion with the previous 
region. We make a mass balance of oxygen concentration G, and we 
assume that all oxygen moles arriving on the inner boundary a(t) react. 
This assumption is not really true, in fact the water plays a key role in 
the speed of reaction, but in our experiments we have a relative humidity 
near to 100%, and this fact justify our assumption; also, the diffusivity D g 
includes implicitly the finite time of reaction. 
The flux of oxygen is 



(25) } g = n p \-D g -^-Gcv p a 



and the mass balance equation is 

dG d 2 G dG 



( 26 ) — - D g-^-Up»— = °- 



The value of G at the boundary j3(£) is given by the value of oxygen on 
the boundary jS, given by eq. |24|) . 

(27) G(j8(f)) = 0(pm 

Since oxygen reacts totally with copper at the boundary a(t), we have 

(28) G(a(t)) = 0; 

Now, we need further condition to close the system (second free boundary). 
Since the oxygen react totally at the boundary a(t), the copper moles wasted 
are given by 



D g dG 1 p. 



c 



(29) 



3.4. Rescaling. The geometry of the problem is given by two regions, one is 
given by x e [/3(f), a(t)] which describes the layer of CM2O, the other is given 
by x € [y(t), /3(f)] i.e. the brochantite layer. Also we have y = - (a> p a + a)},bj. 

3.4.1. New variables. In the first region we adopt (x, t) — > (y, x), in the second 
one (x, t) — > (z, t). We find out 

(30a) y = X — R ) d x = l/(a - p)d y , d xx = l/(a - ^fd XJXJ , d t = f(t)d y + \d % ; 

CI p ty 

(30b) z = % -f^-) d x = l/(j3 - y)d z , d xx = l/(/3 - yfd zz , d t = q(t)d z + \d T ; 
P y t r 
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where 
(31a) 

(31b) 



q(z, t) = 



1 y{d^ - d T a) - d T p 

t r a- $ 

1 z(d T y - d T p) - d T y 



p-y 



Now, we rescale the following parameters and variables in non- 
dimensional form: 



tr 



tr_ 

A 2 ' 



(32b) D g =i^D s , D w = j^D w , D s = -^D s 



From now on, we assume that dotted variables are derivatives with 
respect to the dimensionless time t (V — dV/dx). 

Assuming- O - In if) V - 2--LMhl£l n - An f) Alsn 

Assuming L2 S - ^ n b^s Ms p v ' l w - 2n b M p W r > Li S ~ * n P u gM g p c - Als °' 

y - -((OpU + a)bb), the non-dimensional system is: 
3.4.2. Outer region. For y < x < [} — > < z < 1 



(33) 
(34) 



dS 
dz 



jX d 2 S 

,\ 2 dz 2 



KM 



dS 
dz' 



S(0,T) = S fl 



(35) 
(36) 

(37) 
(38) 



S(1,t) = 0, 



Q s dS 
(hf) dz 

dW _ D w d 2 W _ 
~dz~ ~ To S 2 dz 2 ' 



(1,t) = 6, 



(hff dz2 KM ) 



W(0,t) = W a (z). 



dW 
dz' 



(39) 



10 



F. CLARELLI, B. DE FILIPPO, AND R. NATALINI 



(41) O(0,t) = O„(t). 



(42) ^-)^ (1 ' T) = ^- fc )°-4^a^' 

3.4.3. Inner region. For jS<x<fl— »0<y<l, 

(we have to highlight that fl is a derivative with respect to the t) 
(44) G(l,r) = 0, 



(45) G(0,x)=O 



P' 



(46) X3~ = «- 

From now on we use the non-dimensional system, so we indicate all non- 
dimensional terms without hat. 

3.5. Numerical solutions. To solve our model, we have to set up an 
appropriate numerical scheme which is able to describe the process in 
a short range of time (about 40 hours), but also simulations of 1 year 
taking under control the numerical stability. For these reasons we use 
finite differences schemes with implicit-explicit terms. It can be found in 
II and Il23l . 

3.5.1. Initial conditions. Our numerical procedure requires a(0) + and 
/3(0) + a(0) to avoid singularities in the internal region confining with copper. 
By physical assumption, we know that fl(0) > because it represents the 
first copper consumption. In order that the external equations work (SO2 
is present), we need to set up jS(0) > and y(0) + /3(0). 

3.5.2. Numerical scheme. Indicating explicit term as H(U) and implicit term 
by G(U), we present our system in the external region (eqs. I33l39l) in the 
following form 

(47) U t = H(U) + G(U); 
where 



U = 

the explicit term H(U) is 



( S \ 
W 
O 
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y 



H(U) 



and the implicit term G(U) is 



G(U) 



\dS 
y \ dW 

(p-y) +q j dz 
D s d 2 S 

(fi - yf dzl 

D w d 2 W 

(fi - yf 1* 

D d 2 



, (fi - yf dz 2 

since G(U) is a stiff term which will be integrated implicitly to avoid 
excessively small time steps. 

A more general scheme is IMEX-DIRK Runge-Kutta, which is given by, 
for t - nAt, 

i-1 i 

(48) u® = u n + Af a ik H(u (k) ) + At Y a ik G{u (k) ), i = 1, ...v 



k=l 



(49) 



,n+l 



u n + 



V V 

At ^ &iH{u {i) ) + At Y (OiG{u {i) ). 
i=l i=l 



Here, the matrices A - where % = for / > i and A - (a%) are v x v 
matrices such that the resulting scheme is explicit in H and implicit in G. 
The DIRK formulation requires — for j > i |2|. 

Following the IMEX formalism we shall use the notation Name(s; v; p) to 
identify a scheme, where s is the number of stages of the implicit scheme, v 
is the number of explicit stages and p is the combined order of the scheme. 
In our case we choose to use an Implicit-Explicit Midpoint(l,2,2): s = 1, 
v = 2 and p - 2. 

By the same method, and at the same time-step of previous integration, 
we integrated also equations 



4. Assessment of the model 

4.1. Experiments. The calibration of the model has been made with 
reference to a precise experimental campaign, which is described in 0, 
where the main characteristics of the experimental setting are reported. 
Experiments were carried out in a cyclic corrosion cabinet (Erchisen Mod. 
519/AUTO) and the Cu - \2Sn cast bronze specimens were chosen as a 
representative alloy used in the past (ancient Greek type). For the corrosion 
test, all bronze specimens were exposed to an atmosphere containing about 
200 ppm of SO2 at 40° C and 100% RH for 8 hours (wet cycle), subsequently 
they were exposed to room conditions for 16 hours (dry cycle). Each wet and 
dry cycle was repeated 20 times. The measures of patina thickness, useful 
for the model calibration, were performed using a SEM-EDS instrument 
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Table 1. Parameters table 



Param. 


Value 


Dimension 


Indications 


Pp 


6.00 


g/cm 3 


C112O mass density 


M p 


143.09 


g/mol 


CU2O molar mass 


Pc 


8.94 


g/cm 3 


Copper mass density 


M c 


63.55 


g/mol 


Copper molar mass 


Pb 


3.97 


g/cm 3 


Brochantite mass density 


M b 


452.3 


g/mol 


Brochantite molar mass 


Ps 


1.46 


g/cm 3 


SO2 mass density 


M s 


64.07 


g/mol 


SO2 molar mass 



equipped with an Image Analyzer (IA). Each measurement (after 8, 24 and 
40 hours) has been obtained by an average of 20 measures. 




WD= 15 mm Contrast = 51 .8 % Spot Size = 344 Time :1 6:03:56 



Figure 3. Example of patina thickness measurement with the SEM 
employment. From 171 

The early stage of exposure of bronze samples mainly produced copper 
hydroxyl-sulfate (brochantite and chalcanthite), cuprous oxide, and tin 
sulfide (ottemannite). The tarnish product layer also contained a trace 
of tin oxide. As determined by XRD analysis, the brochantite formation 
was predominant with the increase of the number of corrosion cycles. 
Otherwise, the chalcanthite formation was detected during the early stage 
of the corrosion processes and the intensity of its characteristic peaks 
decreases as a function of time. The same consideration can be made for the 
ottemannite formation, detected with SEM-EDS in localized micro-crack 
areas, which are rapidly covered by basic copper sulfates HZHH). 

4.2. Calibration. The laboratory corrosion tests has been used to calibrate 
the model. We suppose that the parameters to calibrate by experimental 
tests are the diffusivity coefficients. To do that, we used the thickness of 
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corrosion products measured in some time points. Each thickness value has 
been obtained by an average of 20 measures, see figure |2] for two measures 
of them, and we have obtained the following values: 

Table 2. Patina thickness measures. 



Time of measure (hours) 


Averaged value (cm) 


Standard deviation 


8 


5.4418 • 10" 4 


1.7331 • 10" 4 


24 


9.2672 • 10" 4 


1.8473 • 10" 4 


40 


13.2522 • 10" 4 


2.4102 • 10" 4 



Then, we used the least square method to find the best parameters. We 
obtained (in cm 2 /sec): D g = 9.9 • 10" 9 , D s = 3.96 • 10" 5 , D = 9.9 • 10" 6 and 
D w = 3.96 • 10~ 5 . The evolution in time of the corrosion products thickness 
is in figure HI We can see the difference between the experimental points 




Figure 4. Evolution of y(t), /3(f) and a(t). 

and the best simulation in figure |5] 

Remark Here, we test numerical scheme to reproduce the chemical 
reactions in eqs. Q} and ®. Simulations, after 40 hours, give us the 
following values: y = -9.505-10" 4 (cm), a = 3.1693-10" 4 (cm),b = 7.9916-10" 4 
(cm). Knowing the molar volume of copper V c = p c /M c , cuprite V p = p p /M p , 
and brochantite V/, = pb/Mi,, we find that the number of copper moles 
wasted are a/V c — 4.4588 • 10~ 5 , the number of cuprite moles formed are 
hp/Vp = 2.2294 • 10" 5 , i.e. two copper moles are wasted to form one cuprite 
mole, as we expect from eq. Similarly, we obtain the number of cuprite 
moles wasted by reaction |3]l is b/V p - 2.2173 • 10~ 5 , while the number of 
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Experimental points 
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Figure 5. Simulation (continuous line) and experimental points (in red). 

brochantite moles formed by the same reaction is h^/V}, = 1.1086 • 10~ 5 , i.e. 
for each brochantite mole formed two cuprite moles are wasted, as expected 
from ®. 

5. Application 

In this section we assume to use the calibrated model with environmental 
data, detected in Rome at Piazzale Fermi during the year 2005. This 
application is just an example on how the model can be used. Since the 
model has been calibrated for high relative humidity and SO2 concentration, 
the calibration should be improved by experiments with lower values of 
relative humidity and SO2 concentration. The SO2 concentration detected 
is in figure [6] and the temperature is in figure [71 




Figure 6. SO2 concentration detected at Piazzale Fermi (Rome, Italy) 
during the 2005. 

Assuming to have a good calibration with our environmental data, 
we need of saturated vapor density (SVD in [g/cm 3 ]) as function of 
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Figure 7. Temperature detected at Piazzale Fermi (Rome, Italy) during 
the 2005. 

relative humidity and temperature T [°C] detected at Piazzale Fermi. It 
is useful for getting an exact quantity of water vapor in the air from 
a relative humidity (RH), the density of water in the air is given by 
RH ■ SVD - ActualVaporDensity. The behavior of water vapor density 
is a non-linear function, but an approximate calculation of saturated vapor 
density can be made from an empirical fit of the vapor density curve [50] 
(between - 45 °C is a good approximation). 

(50) SVD(T) = 5.018 + 0.32321T + 8.1847 • 10 _3 T 2 + 3.1243 • 10 _4 T 3 . 
The figure of saturated vapor density obtained is in figure |H] 




Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 



Figure 8 . Saturated vapor density obtained by temperature and relative 
humidity detected at Piazzale Fermi (Rome, Italy). 

Thus, using SO2, temperature and humidity data detected at Piazzale 
Fermi in Rome during the 2005, we obtained the behavior of corrosion 
products illustrated in figure |H 
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Figure 9. Simulation of corrosion with environmental data detected at 
Piazzale Fermi (Rome, Italy) in 12 months (2005). 

Results of this simulation represent a qualitative example of corrosion 
formation on copper under SO2 attack. They have to be improved by others 
laboratory experiments, because our model is calibrated using a very high 
SO2 concentration (never present in the atmospheric environment in this 
concentration), and under a very high relative humidity (near to 100%). 
To have a more accurate calibration, it would be necessary set up mainly 
experiments with different SO2 concentration and different concentration 
of humidity in the air, so to calibrate the model under these variations. 

Our conditions produce a very thick cuprite layer in 1-year simulations 
with respect to the laboratory experiments, because of the air SO2 
concentration is very lower than that present in the experimental room, 
and this caused a slower growth of brochantite layer. 
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6. Conclusions 

In conclusion a new mathematical model was developed to describe 
and simulate the evolution of brochantite formation on Cultural Heritage 
artifacts exposed to sulfur dioxide atmosphere. The aim was to create a 
new approach to forecasting corrosion behavior without the necessity of an 
extensive use of laboratory testing using chemical-physical technologies, 
while taking into account the main chemical reactions. Although the model 
was kept simple, just describing the main reaction and transport processes 
involved, the mathematical simulations and the related model calibration 
are in agreement with the laboratory experiments. 
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